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Abstract 


This  report  treats  the  description  of  atomic  systems  by  the  novel 
method  of  introducing  interelectron  coordinates  explicitly.   The  simplest 
atomic  systems,  heliijm  and  lithium,  are  treated  and  excellent  results  are 
found  in  calculating  the  energy  values  of  the  corresponding  ground  states. 
In  addition,  an  approximation  technique  of  high  accuracy  has  been  developed 
to  handle,  mathematically,  wave  functions  containing  many  interelectron 
coordinates . 
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Introduction 


The  problem  of  suitably  describing  atomic  systems  has  been  attacked 
since  the  advent  of  Quantum  Theory.  The  statistical  theory  of  Thomas'-  -^ 
and  Fermi  ^  -' ^  the  self-consistent  method  of  Hartree '-  -'  and  the  methods 

M  [5! 

of  Slater '-  -'  and  Morse"-  -"  are  among  the  oldest  and  have  been  widely 
investigated.   All  of  these  methods  deal,  in  essence,  with  the  nuclear 
coordinates  of  the  individual  electrons.   Although  it  has  long  been 
recognized  that  the  interelectron  potential  within  the  atom  is  of  com- 
parable magnitude  to  the  nuclear  potential,  little  work  was  done  original- 
ly in  using  interelectron  coordinates  explicitly.   The  well-known  varia- 
tional method  of  Hylleraas  ^  -I  does  introduce  the  interelectron  coordinates. 
However,  the  wave  function  is  specified  only  after  a  variational  calcula- 
tion of  the  energy,  and  the  method  is  consequently  difficult  to  apply  to 
excited  states  and  to  atoms  much  heavier  than  lithium.   The  method  of 
configuration  interaction  as  generalized  by  Lowdin '-  -'  and  applied  by 
Tycko,  Thomas  and  King'-  -'  shows  excellent  promise  in  treating  the  inter- 
electron effects  although  the  interelectron  coordinates  are  not  intro- 
duced  in  a  direct  manner.   The  recent  approaches  of  Brueckner'-  -■  and  of 
Bohm  and  Pines'-  -^  which  also  treat  interelectron  effects  have  not  found 
much  application^  as  yet,  in  atomic  systems. 

A  few  years  ago,  Pluvinage '-  -'  made  a  substantial  advance  in  the  use 
of  interelectron  coordinates  in  atomic  systems  when  he  illustrated  how  the 
Schroedinger  Equation  can  be  partially  separated,  using  nuclear  and 
interelectron  coordinates  together,  so  that  the  interelectron  potential 
no  longer  appears  as  the  perturbing  term.   In  this  partial  separation 


the  perturbing  term  has  the  character  of  derivative  terms  (as  it  were^ 
'velocity  dependent'  potentials)  and  these  derivative  terms  are  always 
finite  in  space.   The  interelectron  potentials  are,  of  course,  infinite 
when  two  electrons  have  zero  separation.   The  zero  order  wave  function 
for  such  a  system  becomes  a  product  of  bound  state  hydrogenic  wave 
functions  for  the  nuclear  coordinates  and  continuum  state  hydrogenic 
wave  functions  of  variable  kinetic  energy  for  the  interelectron  coordin- 
ates.  Such  a  wave  function  takes  into  account  the  binding  of  all  the 
electrons  to  the  nucleus  and  also  the  mutual  scattering  among  the  electrons 
themselves.   We  will  designate  product  wave  fimctions  of  this  type  as 
scattering  wave  functions  (SWF)  to  distinguish  them  from  the  product  of 
bound  wave  functions  (BWF)  involving  only  the  nuclear  coordinates. 

The  SWF  would  be  expected  to  yield- much  better  results  than  the 
corresponding  BWF  since  all  of  the  potentials,  nuclear  and  interelectronic, 
are  accounted  for  in  the  SWF.   This  was  indeed  the  result  found  by 
Pluvinage  when  he  used  the  SWF  to  obtain  variational  energies  of  ground 
states  of  helium- like  systems.   The  SWF  gave  answers  which  were  in 
general  twice  as  close  to  the  experimental  results  as  the  corresponding 
BWF. 

The  purpose  of  this  paper  is  to  modify  and  expand  the  Pluvinage 
approach  in  such  a  way  that  SWF  can  be  used  with  a  reasonable  amount  of 
labor  for  atomic  systems  beyond  the  helium  type.   Once  the  modification 
has  been  accomplished  the  actual  results  will  be  tested  in  the  case  of 
the  ground  state  of  lithii:m. 


Two  things  have  been  done  in  applying  SWF  to  atomic  systems  beyond 
helium.   It  vas  first  noted  that  the  use  of  continuum  vave  functions  of 
zero  kinetic  energy  greatly  simplifies  the  mathematics  ajid  allows  integra- 
tion of  the  SWF  terms  of  tabulated  functions.   In  Section  2  we  discuss  this 
aspect  of  the  problem  and  present  a  comparison  of  the  energies  calculated 
using  the  simplified  SWF  with  those  of  Pluvinage  in  helium-like  systems. 
In  Section  3  the  method  for  applying  these  simplified  SWF  to  more  compli- 
cated systems  is  described  and  it  is  seen  that  these  SWF  describe  the 
nuclear  coordinates  as  Fermi-Dirac  'particles'  filling  the  bound  energy 
levels  pairwlse,  while  the  interelectron  coordinates  are  described  as 
Bose  'particles'  all  piling  into  the  lowest  energy  continuum  state. 

The  second  thing  done  in  applying  SWF  was  to  devise  an  accurate  approx- 
imation method  to  integrate  functions  containing  three  or  more  interelectron 
coordinates  over  the  space  of  the  nuclear  coordinates  (including  angles). 
In  Section  h-   an  accurate  approximation  method  is  described  which  consists 
of  expanding  the  functions  containing  the  interelectron  coordinates  about 
well-chosen  points  in  the  space  of  the  nuclear  coordinates.   The  approx- 
imation technique  is  tested  in  helium-like  systems  and  is  found  to  be 
better  than  99  /o     accurate. 

In  Section  5,  the  ground  state  energy  for  lithium  is  calculated  using 
the  SWF  proposed  and  the  approximation  technique  described.   The  accuracy 
of  the  approximation  technique  is  fully  maintained  in  lithium.   The 
results  are  discussed  to  some  degree  in  Section  6. 


k  - 


2.    Scattering  vave  fimctions  in  helium-like  systems 

Pluvinage  treated  the  problems  of  two  electrons  in  the  ground  state 
of  helium-like  systems.   The  wave  equation  for  such  a  two  electron  system 
is  the  following: 


(2.1)  H\|/  = 


W-%v.2)    -^-^.1 


2      1  2  r  r  r 
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ilf  =  E\lf 


using  atomic  units.   The  notation  has  its  customary  meaning.   Because  of 
the  Coulomb  potentials,  the  wave  equation  possesses  singularities  when 
the  electronic  distances,  r  ,  r  ,    are  zero  or  when  the  interelectron  spac- 
ing, r,  ,  is  zero.   The  unique  feature  of  Pluvinage 's  approach  is  that  he 
separated  the  Hamiltonian  in  such  a  way  that  the  behavior  of  the  electrons 
at  any  of  the  singularities  is  exactly  correct.   Consequently  the  effect 
of  all  the  potentials,  nuclear  and  interelectron,  is  taken  into  account 
to  a  very  large  extent. 

If  the  coordinates  r  ,  r  ,    r   are  used  and  symmetric  states  only 
are  considered,  then  H  becomes 
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(2.5)     2  cos  0. 
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and  9.  . .  Is  the  angle  between  r.  and  r. .  . 

t 
This  result  inmediately  suggests  that  we  take  H  =  H  +  H  with  the 

unperturbed  wave  equation  as 
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The  solution  is  the  SWF 


(2.5)      ^„  =  ^^ir^)^J-2K^^l2^ 
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Here  0  is  the  hydrogenic  wave  function  for  the  n   bound  state  of  an  atom 


whos 


e  nuclear  charge  is  Z,    \r     is  the  hydrogenic  wave  function  for  an  electron 

o 
in  the  continuum  having  an  energy  k  /2  and  scattering  from  a  center  of 

charge,  l/2.   The  factors ;,  l/2,  occur  because  the  interelectron  coordinate, 

r  _,  is  shared  by  two  electrons.   The  perturbation  is  now 


(2.6)     H  =  -  cos  0 
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H  is  finite  for  all  values  of  the  coordinates  and  is  small  compared  to  H 


at  the  singularities.   The  solution  to  the  unperturbed  problem  thus  repre- 
sents the  actual  wave  function  exactly  at  all  the  singularities. 
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Pluvinage  took  the  SWF  from  Equation  (2.4)  and  used  it  to  calculate 
the  energy  of  the  ground  state  by  substituting  it  in  the  variational 
expression 


(2.7)  E  =  -^ ^  =   -  Z^+   k2+  -^ ^ 


The  integrals  were  evaluated  "by  power  series  and  the  parameter  varied  was  k. 
This  resulted  in  the  energy  value^  -2.873  AU,  for  helium.   The  correct 
energy  is  -2.90572  AU,  while  the  energy  obtained  by  using  the  product  of 
bound  hydrogenic  wave  function  with  variable  atomic  number,  Z,  is  -2.848  AU. 
The  SWF  function  yields  an  evergy  value  which  is  closer  to  the  true  energy 
by  a  factor  of  two  when  compared  to  comparable  BWF.   This  result  is  also 
true  for  the  calculations  made  upon  other  helium-like  systems. 

The  good  energy  value  calculated  by  the  Pluvinage  approach  is  grati- 
fying. Even  if  no  better  energy  value  were  found  the  approach  deserves  a 
great  deal  of  consideration  as  a  fundamentally  different  approach  to  atomic 
systems.   The  treatment  of  the  singularities  due  to  the  potentials  in  the 
wave  equation  is  very  satisfying  and  lends  a  great  deal  of  validity  to  the 
method. 

The  true  wave  function  for  the  lowest  state  of  a  two  electron  system 

is  expected  to  have  nodes  only  when  the  electrons  are  at  infinity  and  at  the 

nucleus.   However,  the  wave  function  (2.5)  has  an  infinite  number  of  nodes 

when  k  >  0  since  u,  oscillates  at  large  distances.   On  the  other  hand,  u 

k  o 

has  no  nodes  at  any  value  of  the  argtmient.   In  addition,  physical  grounds 
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would  indicate  that  k  should  be  zero.   If  we  removed  the  two  electrons  far 
from  the  nucleus,  their  kinetic  energies  would  decrease  and  their  motion 
relative  to  each  other  would  also  decrease  and  at  an  infinite  distance 
from  the  nucleus  they  would  scatter  from  each  other  with  no  relative 
motion  at  all. 

The  unperturbed  wave  equation  for  the  interelectron  coordinate  with 
k  =  0  is 


2 
d  u      2du       u 

(2.8)  1  +  -^ 2 


dr        dr       r 
12        12      12 


The  solution  of  this  equation  is  the  modified  Bessel-Clifford  function  of 
the  first  order,  E  ;,  which  is  finite  at  r  ^  =   0.   Thus 


(2.9)  u^  =  — —:  J^  (21 J^)^     E^(r^) 

^12 


V- 


where  J  is  the  Bessel  function  of  first  order.   Tables  of  the  Bessel- 


Clifford  function  are  available 


[12] 


The  advantage  of  the  use  of  the  wave  function,  u  =  E  ,  for  the  SWF 
is  twofold.   Firstly,  the  integrals  which  contain  the  full  wave  function 
are  integrable  in  closed  and  relatively  simple  form.   Secondly,  a  physical 
basis  has  been  obtained  for  selecting  one  of  a  set  of  SWF  for  use  in  pro- 
blems.  This  obviates  the  need  to  adjust  k  by  a  variational  method  when- 
ever the  SWF  is  to  be  used.   This  second  is  a  true  advantage  only  if  the 
wave  fimction  using  k  =  0  is  comparable  in  accuracy  with  one  using  the  best 
value  of  k.   That  this  is  so  is  shown  below.   The  a  priori  selection  of  a 
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particular  value  of  k  will  save  a  considerable  amount  of  effort  in  the 
calculation  of  the  energy  level  of  the  lithium  ground  state. 

We  now  proceed  to  the  calculation  of  the  energy  of  the  ground  state 
of  helium  and  helium-like  ions  using 


-Zr   -Zr 
(2.10)         M;^  =  e   ^  e     ^^(^12^ 


This  will  be  compared  with  Pluvinage ' s  results  as  a  test  of  the  useful- 
ness of  this  particular  SWF.  As  usual  the  variables  used  by  Hylleraas '-  -' 
are  best  suited  to  the  problem.   Therefore^,  we  set 


(2.11)         s  =  r^+  r^  ,    t  =  r^-  r^ 


and  retain  the  variable  r   •   The  volume  element  is  now 


(2.12)         dt  =  2rt^(s^-  t^)r  ds  dt  dr 


Since  the  ground  state  is  symmetric,  only  even  powers  of  t  appear  in  the 
integrals  and  this  results  in  the  factor  2  appearing  in  the  equation.   It 
is  very  convenient  to  arrange  the  limits  of  integration  in  such  a  manner 
that  the  integration  over  r   appears  last  and  covers  the  region  from 
zero  to  infinity.   Thus  we  have 


(2.13)         0  <  t  <  r^2  ,    1-12  -  ^  -  °°  ^    0  <  ^12  5  °° 
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Using  the  Hylleraas  variables  and  remembering  that  k  =  0^  the  ex- 
pression for  the  energy  becomes 


(2.1I1)  E  =  -  Z^  +  ^ 

00 


vhere 


00      00    r^2  /  dE  ^ 

(2.15)    H^:  =  Z(2.2)  j     a^^^j     as  /   dt  s(r^|-t2)  e-2^  f^^ 
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Since  k  is  constant  the  variational  expression  for  the  energy  has  become 

the  usual  first  order  perturbation  expression. 

The  quantities  W   and  H   are  evaluated  by  first  integrating  over 
^  00      00  '' 

t 
the  t  and  s  coordinates.   In  the  case  of  H   an  integration  by  parts  is 

then  performed  and  ve  have  finally 


(2.17)    H„;  =  2„2z(|l,(Z)-il3(Z)-^yz)) 

(a.i8)        N^^  =  2„^  ig  (  f  i,(z)  .  I  Ij(z)  .  -ij  I^(Z)) 


10  - 


vhere 
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:ionM 
with  the  substitution  r   =  V^.   The  expression  for  I,  is 


The  integral^  I  ,    is  evaluated  by  integrating  over  the  Bessel  functi 


(2.20)     I^(Z)=  ^exp(|)E^(-i5) 


itz'    ^'"    "Hz' 


The  integrals  for  n  greater  than  one  are  obtained  easily  from  I  by  differ- 
entiating with  respect  to  Z  or  by  using  certain  recurrence  relations.   The 
method  for  doing  this  is  discussed  in  Appendix  I  where  the  explicit  expres- 
sions for  several  of  the  I   are  also  given.   These  expressions  are  simple 
combinations  of  polynomials  in  Z,  the  exponential  function  and  the  modified 

Clifford- Bessel  functions,  E  and  E,  .   The  latter  are  rather  fully  tabulated 

o      1 

in  Reference  [12] .   They  can  also  be  computed  from  their  power  series 

representation,  Equation  (1.3)^  which  converges  quite  rapidly. 

The  ground  state  energy  was  calculated  for  several  helium-like  ions 

and  atoms  by  the  method  outlined  above.  Its  deviation,  D  ,  from  the  experi- 

mental  energy,  (D  =  E     —  E    )  is  given  in  Table  I,  in  atomic  units,  and 

calc    exp 

is  compared  with  several  other  computed  deviations.  D  is  the  deviation 
obtained  from  Pluvinage's  work,  D  is  that  calculated  using  BWE  and,  D 
that  calculated  by  BWF  with  variable  Z. 
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There  are  several  noteworthy  features  exhibited  in  Table  I.   The 
indicated  superiority  of  the  Pluvinage  wave  function  is  clearly  shown.   That 
the  choice  of  k  equal  to  zero  is  indeed  a  good  one  is  quite  evident.   A  major 
part  of  the  improvement  in  energy  obtained  by  Pluvinage  is  retained  with  k 
equal  to  zero  with  a  substantial  saving  in  labor.   The  energy  value  of  H  is 
easily  calculated  with  k  taken  as  zero  whereas  Pluvinage  did  not  calculate 
this  energy,  presumably  because  of  the  slow  convergence  in  the  power  series 
he  used.   Incidentally,  the  energy  value  calculated  for  H  does  not  quite 
yield  binding,  falling  0.002  AU  short  of  the  energy  of  atomic  hydrogen. 

Since  the  two  electrons  in  the  inner  shell  of  neutral  lithium  contri- 
bute the  major  part  of  the  energy  in  the  ground  state  of  lithium  (  Z  =  5), 
the  table  shows  that  the  use  of  E   in  the  wave  function  will  probably  yield 
80  %  of  the  improvement  we  can  expect  if  we  varied  k.   The  labor  saved  in 
the  process  of  calculating  the  energy  of  neutral  lithium  using  E  will  prove 
substantial. 


D  =  E   ,   -  E 

calc    exp 

in  A.U. 

z 

E    (A.U.) 
exp^ 

^B 

\ 

^S 

^p 

h" 

1 

-0.528 

0.155 

0.055 

0.050 

_ 

He 

2 

-2.90i+ 

0.154 

0.056 

0.029 

0.026 

Li+ 

3 

-7.280 

0.155 

0.057 

0.051 

0.025 

Be2+ 

h 

-15.656 

0.155 

0.057 

0.055 

0.025 

b3+ 

5 

-22.051 

0.156 

0.058 

0.054 

0.025 

ch+ 

6 

-52.406 

0.156 

0.058 

0.054 

0.025 

n5+ 

7 

-44.781 

0.156 

0.058 

0.055 

0.025 

00 

- 

0.157 

0.059 

o.o4o 

~ 

TABLE  I. 


Deviation  of  various  calculated  energy  values  from  the  experi- 
mental values  of  the  ground  state  of  helium-like  atoms  and  ions. 
The  energy  differences  at  infinite  Z  are  obtained  by  comparison 
with  the  formula  given  by  Hylleraas  and  Mitdal  [l4j  for  the  energy. 
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The   advantage   of  using  the  SWF  vith  k  =  0,    over  the  BWF  in  a  per- 
turbation-type  approach  to  atomic  wave   functions  is  even  more  evident  in 
Table  I.      D  ^  is  the  deviation  of  the   first   order  perturbation  energy 
from  the  experimental  value  vhen  the   interelectron  distance   is   considered 
the  perturbation.      This   deviation  is   four  to  five  times  the   deviation 
calculated  vhen  the  Hamiltonian  is  broken  in  such  a  manner  that  Equation 
(2.6)   is  the  pertiirbation.      This  strongly  Indicates  that  excited  states, 
vhlch  cannot  be  treated  easily  by  a  variational  method,    could  be  treated 
■with  good  results  by  a  perturbation  approach  based  upon  the  SWF  vith  k 
taken  as   zero. 


5.       Many-electron  systems 

In  the   last  section  ve   illustrated,    in  the   case  of  the  tvo- electron 
system,    the  simplification  ve  have  proposed  to  the  Pluvinage  method.      The 
extension  to  the  many  electron  case  poses  some  problems.     We   cannot  vrite 
dovn  a  Hamiltonian  for  a  many  electron  system  vhlch  is  analogous  to  the 
one  that  has  been  vrltten  down  for  the   case   of  helium.      The  easiest  vay 
to  see  this  is  to  note   that   if  ve  wei^  to  introduce  the   interelectron  co- 
ordinates,   the  radial  coordinates   of  each  electron  plus  three  Euler  angles, 
the  number  of  coordinates  we  have   available  with  which  to  express  the 
Hamiltonian  is 

N  +  II(N  -  l)/2  +  3 

where  N  is  the  number  of  electrons.  The  original  number  of  degrees  of 
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of  freedom  is  JW.   If  we  nov  insist  that  the  total  numher  of  coordinates 
used  to  express  the  Hamlltonian  be  the  same^  we  have  that  N  must  be  equal 
to  either  2  or  5.   Thus,  ve  could  not  proceed  in  this  vay  for  any  element 
beyond  lithium. 

We  can  circumvent  this  difficulty  as  follows:  First  we  are  guided 
by  our  discussion  of  the  helium  problem  in  our  choice  of  an  approximate 
wuve  function  to  describe  the  many  electron  system.  The  generalization 
of  the  wave  function  which  we  have  used  for  the  helium  atom  is 


(5-1)  ^o 


i        J  /  >  n 


J]     ^,ir.) 


in  vhich  D  represents  the  Slater  determinant  of  unperturbed  hydrogenic 
vave  functions   corresponding  to  suitable  unperturbed  states.      The  E 
functions  are  the   Cliff ord-Bessel  functions   of  their  argument.      The  pro- 
duct of  the   Cliff ord-Bessel  function  remains   invariant   as  we   interchange 
the   coordinates.      Thus,    the   total  vave  function  is  properly  antisymmetric 
provided  we  include  the   spin  wave  functions   in  the   0    .      According  to   (5«l) 
we  have  restricted  ourselves   to  S   states  since  the   0   ^s  depend  only  on 
the  radial  coordinates.      The  generalization  to  other  unperturbed  angular 
momentum  states  poses  no  essential  difficulty. 
The  approximate  energy  is  given  by 

(5.2)  (    \lf      ,    E  ^      )      =      1    \lr*  H  t     drT"  .    .    .    d^ 

^  ^      o  o  tool  n 


ll^  - 


where  H  is   given  by 


(3.3)  \-l^\' 


\|f    =   E\(r 


and  the   \(f     are  properly  normalized.      Since  the  \|f     depends  on  the   coordinates 

r,    and  r. ..   we   can  reexpress  the   operator  H  in  terms  of  the   coordinates  r. 

i  ij  i 

and  r       by  the   rules  of  partial  differentiation.      This  deccmposes  H  into 
two  parts  H     and  H   ,   where 


(3.M  H^ 


1  /   _J_        d  2d 

2  I        2     "S??     ^i       "SrT 
^i  ^  ^ 


r„       -s + 


/>m     \      7m  ^^  -^^  ^^ 


and 


(3-5)  h'   =  -    XI      H     cos  0.'^. 


a^ 


Sr . Sr .  . 


H  E  E 


3= 

COS   9, 


a?he  effect  of  H  operating  on  \lf  in  the  integral  is  to  yield  a  contribution 

to  the  energy  equal  to  the  sums  of  the  hydrogenic  energies  of  the  electrons 

neglecting  their  interactions,  since  the  use  of  the  Cliff ord-Bessel  functions 

in  \lf  corresponds  to  an  eigenvalue  of  0  for  that  part  of  H  • 
o  o 
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The  corrections  to  the  hydrogenic  energies  is  given  "by 


*   ' 


(5-6)  (  t^  H  1^^  )dr^  .  .  .  dr^ 

The  perturbing  Hamiltoniaxi,  H  ,  contains  a  nev  set  of  terms  vhich 
do  not  appear  in  tvo  electron  systemB.  In  addition  to  the  interaction 
between  the  nuclear  and  the  interelectron  coordinates,  the  interaction 
among  the  interelectron  coordinates  now  appears. 

When  we  had  to  evaluate  the  integral  similar  to  (5.6)  for  the  case 
of  helium,  we  expressed  the  volume  element  in  terms  of  the  coordinates, 
r^^,  r^  and  r^,  and  properly  restricted  the  range  of  integration.  For 
heavier  atoms  expressing  the  volume  element  In  terms  of  the  interelectron 
coordinates  is  not  a  fruitful  way  to  proceed.   In  the  first  place,  there 
axe  too  many  of  them  and,  in  the  second  place,  it  is  doubtful  that  the 
integrals  could  be  carried  out  since  the  range  of  integration  would  have 
to  be  severly  restricted.  However,  this  paper  will  present  a  very 
accurate  approximation  technique  for  doing  this  in  a  rather  simple  manner 
and  this  technique  will  form  the  basis  of  the  calculations  in  lithium. 
The  lithium  problem  will  thus  serve  as  a  testing  ground  for  the  useful- 
ness and  accuracy  of  the  approximation  technique  developed. 
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k.        Mathematical  approximation 

The  type  of  integration  which  must  be  evaluated  is  of  the  follow- 
ing general  foim: 


(^•1)     I^^^  =/  ^^^'    ^1'    ^2 ^n)s^^l2'  ^15----Vl,n) 


dV 


The  function,  f,  has  a  monotonic  exponential  variation  in  the  radial 
direction,   g  also  varies  monotonically  but  more  slowly  than  f.  dV  is  the 
volume  element  in  the  space  of  the  n  electrons.   The  integral  of  f  alone 
over  all  the  coordinates  can  be  performed  and,  in  the  applications  we 
make,  is  simply  the  integral  of  products  of  hydrogenic  wave  functions 
over  the  corresponding  radial  coordinates  and  angles.   The  relatively 
slow  variation  of  g  compared  to  f  leads  us  to  expect  that  a  Taylor  expan- 
sion of  g  about  well-chosen  points  in  the  space  of  all  the  radial  coordin- 
ates and  angles  will  lead  to  a  rapidly  convergent  series  of  Integrals  to 
represent  I.   Indeed  we  have  been  able  to  treat  the  expansion  in  such  a 
way  that  the  first  term  alone  gives  an  answer  which  is  generally  more  than 
99  /»  accurate  for  the  type  of  integrals  with  which  we  must  deal.   Such 
accuracy  will  give  an  answer  for  the  energy  of  the  lithium  ground  state 
with  a  computational  error  much  smaller  than  the  expected  difference 
between  the  computed  energy  and  the  experimental  energy. 

The  scheme  adopted  to  evaluate  I  is  as  follows:   Expand  g  about  the 

radial  points  r, „,  r^-  ...  r   thus 
^  10'   20      no 


(if. 2)     g(r.,,  r.^...r_-   )  =  g  +  51  — 


12'  ^15"-  n-l,n'  -  ^o  "  ^—  dr. 

1=1    ^ 


(r.  -  r,  )  + 

^1    10 


r .  =  r . 
1   io 
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vhere  g  indicates  g  evaluated  at  the  radial  expansion  points.   Since  we 


have  not  expanded  g  in  terms  of  the  angles^  g  and  - — 

o     rtr. 

1 


are,  as  yet. 


r.  =  r. 

1     lO 


fionctions  of  the  various  angular  coordinates.   Substitute  Equation  (^.2) 
into  (^.l);  then  there  results 


n 


(If. 5)     I(Z)  =  J  ...  J  fg^dV  +  ^  I  f 

i=l 


dg 


dr. 
1 


(r.-  r.  )dV  + 

1     lO 

r.=  r. 

1    lO 


To  speed  the  convergence  of  the  series  the  r.   are  chosen  such  that  the 
•^  lO 

second  terms  in  the  series  are  zero.   Thus  the  formulas  for  the  r.   are 

00 

(k.k)  \      f  X   (r^-  r^^)dY   =  0  ;   i  =  1,  . . .  n  . 


The  separation  of  the  angular  and  radial  coordinates  alloved 


dr. 

1 


to  be 


taken  outside  of  the  integral. 


r. 

lO 


After  the  r.   have  been  determined  the  angular  expansion  points 
are  then  found.   The  angular  treatment  is  different  from  the  radial  treat- 
ment for  two  reasons.   First,  the  angular  variation,  which  appears  in  the 
term  fdV,  is  generally  slow  or  only  comparable  to  the  variation  that 

occurs  in  g.   Second,  the  radial  treatment  depends  only  on  the  value  of  g 

de 
and  not  its  shape,  since  t-^—  does  not  appear  explicitly  in  determining 

or. 

1 

r.   or  I.  We  make  use  of  the  angular  treatment  to  introduce  the  shape 

10 

of  g.   The  criterion  established  to  determine  the  angular  coordinates 

0.  ,a.   is  then  this:   The  angles  are  chosen  to  be  those  angles  for  which 
10   lo  '^  ^ 

the  integrand  has  a  maximum  subject  to  Equations  (4.ij-). 
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It  is  more  convenient  to  deal  with  the  logarithm  of  the  integrand; 
therefore^ 


(^.5; 


dQ. 


^  v.(0.)g  I 


=   0 


dQ. 


^"  "of«j'«o 


w  is  the  appropriate  weighting  function  for  the  angle  which  appears   in 
the  volume  element. 

The  system  of  Equations  (^.5)  must  be  solved  simultaneously,  and 
in  practice  an  iteration  procedure  is  used  to  solve  them.   Since  the  inter- 
electron  coordinates  which  enter  into  g  do  not  depend  strongly  upon  the 
angles,  the  iteration  of  Equations  (^.5)  converges  quite  rapidly  if  a 
reasonable  guess  is  originally  made  for  the  angles. 

Writing  r. .   as  the  value  of  r. .  at  r.  ,  0.  ,  a.   we  now  have  for  I 
ijo  ij     10   10   10 


{k.6)  I^(Z)  =  g(r^20'  ""150'  ••■    '     Vl^no^  J  ^^^'^1'  ^2^  •*•  '    ""-^    ^^ 


n 


Note  the  dependence  of  I  upon  Z  in  Equations  (^.l)  and  (4.6).  With  large  Z 
the  exponential  variation  of  f  is  so  strong  that  only  those  regions  near 
the  nucleus  contribute  to  I.   As  long  as  g  is  monotonic  this  means  that 
Equation  (4.1)  reduces  to 


(^.7) 


1(00)  =  g(o,o, o)  I  ..  I  f  dV  . 
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This  is  exactly  the  form  vhich  Equation  (^.6)  assumes  at  large  Z  and  that 
expression  is  thus  exact  at  Z  equal  to  infinity. 

The  approximate  expression  for  I  can  be  modified  in  several  ways  of 
special  importance.   For  example^  it  sometimes  happens  that  g  is  not 
monotonic  or  uniquely  defined  when  all  interelectron  coordinates  are  zero. 
However,  it  may  be  possible  to  separate  g  into  two  parts,  one  of  which  is 
monotonic  and  uniquely  defined  at  zero  interelectron  distances  and  the 
remaining  function  is  directly  integrable.   Thus 

(i|.8)  g  =  Gy 

where   y.  is   integrable.      Then,    in  the   same  manner  as   previously. 


(^•9)  I^Z)    =   G(r^2Q,  ....)    I    fxdV 


still  subject  to  Equation  (^.U).  This  expression  is  exact  at  Z  =  oo.  Other 
modifications  are  also  possible.  For  instance,  %  may  be  integrable  when  the 
r.   are  determined,  or  perhaps  it  is  not  integrable  at  all.   In  these  two 

lO 

cases  equations  paralleling  (^.6)  and  (^.9)  are  easily  written;   however, 
they  are  not  correct  at  Z  =  oo .   To  make  them  correct,  the  equations  can 
be  multiplied  by  a  constant  which  equals  the  ratio  of  the  approximate  to 
the  true  answer  as  Z  approaches  infinity. 

The  accuracy  of  the  method  described  above  will  be  illustrated  by 

using  it  to  calculate  the  normalization  integral  N   and  the  perturbing 

t 

Hamiltonian  integral,  H   ,  for  helium- like  systems  and  comparing  with 

the  exact  results  given  by  Equations  (2.l6)  and  (2.17). 
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The  normalization  Integral  is: 


(   r     -2Z(r  +r  )  _^   _^ 

(4.10)        \o=JJ^  ^1  ^^12)  ^^1  ^^2 

where 

-2Z(r^+r2)         ^ 


f  -  e 


,   g  =  E^"(r^2) 


Take   r     as   a  polar  axis   and  let  0     be  the   azimuthal  angle  between  r     and  r   . 


The 


Integrals   over  all  the   orientations  of  r     and  over  the  polar  angle,    0   , 


or  r     can  be  performed  immediately.      This  results   In 

00   001^     „„/  \ 

o    (     f    C    -2Z(r  +r    )  2  2 

(4.11)      N^^  =  8«2  J     j   j  e  ^     ^     \  (r^)    r^  dr^  r^  dr^  sin  9^  dG^ 


000 
where  now 


^12   =    J^l^  4-   2^1  ^2   "°"   ®2 


The   radial  expansion  points   are  determined  by  Equation   (4.4)   yielding 

The   angular  expansion  points   reduce  to  the   single  point,    0p^,    determined 
by  Equation   (4.5)    in  the  form 

(4.15)  ^   [in  (sin  02  E^^  (^^^,    ^^^^    g^)j  ^       . 


Set 
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(^.IK 


0° 


r  dE 

120  1 


dr 


12 


■120 


"o   ^^120^ 
h   ^^120^ 


-   1 


where  we  have  used  Equation   (l.^)    of  the  Appendix.      Thus 


(^.15) 


cos   9^=- 


2  +  0' 


-   1 


and 


(^.16) 


r  =    v/2   r  /2   -   E    /E 

120      ^      10  >y         re 


where  E  and  E^  are  evaluated  at  r^^-. 
o      1  120 


Using  this  value  of  r    and  inte- 


grating over  the  exponential,  the  formula  (^.6)  becomes 


(It. 17) 


C'("'/2')\'(>-X20' 


The  evaluation  of  r    from  (k.l6)    is  obtained  by  a  series  of  itera- 
tions, which  as  mentioned  already,  converges  rapidly.   Two  iterations 


generally  yield  r    to  an  accuracy  of  a  part  per  thousand  with  a  fair 


estimate  for  9   or  r    initially. 
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2 
The  quantity  E   is  a  direct  measure  of  the  change  produced  in  the 

normalization  integral  because  of  the  presence  of  the  interelectron  term 

in  the  wave  function.   A  comparison  of  the  approximate  ansver  for  the 

2 
normalization  integral  with  the  exact  is  made  in  Table  II  where  E 

Z  N   /rt   is  tabulated  versus  Z  N   /rt   for  various  values  of  Z.   Except 
00  '  oo  ' 

for  Z  =  1,  the  error  is  less  than  one  percent  despite  the  fairly  large 
change  produced  by  the  presence  of  the  interelectron  term.   The  inter- 
electron term  is  clearly  not  a  small  perturbation  even  at  Z  as  large  as 
6.   The  excellent  accuracy  throughout  the  entire  range  is  quite  gratify- 
ing;  at  Z  =  1,  the  error  is  less  than  6  /o   despite  the  Increase  of  almost 
9  times  in  the  normalization  integral. 

The  integral^  H   ,  is  obtained  from  (2.6)  by  integration.   It  is 
convenient  for  the  present  case  to  differentiate  with  respect  to  the 
radial  coordinates  and  then  to  combine  the  angular  terms  by  using  Equa- 
tion (2.3).   Establishing  r  as  the  polar  axis  and  integrating  over  the 
orientations  of  r  and  the  polar  angle  of  r  as  before,  yields: 


(i^.i8)    H    =  (at^)z 

00        ^ 


00    00    It 

\\\ 

000 


■2Z(r^+r2) 


(r^+   r^)    (1   -    cos   Qg) 


1       dE 


X 


^12     ^^12 


12  2 

—     r^     dr^  r^     dr^  sin  0^  dO^ 
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The  function  f  is  the  same  as  before  but  e  has  the  form 


(^.19) 


(r^+  v^){l   -    cos  0^)   du^ 


"12 


dr 


12 


Because  of  the  angular  term,  g  is  not  uniquely  defined  at  r   =  0.   We 
thus  try  Equation  (^.8)  where 


dE. 


(4.20) 


G  = 


X 


dr 


12 


(r^+  r2)(l  -  cos  O^] 


"12 


and  we  find  that  %  is  indeed  integrable. 

As  before  r   =  r   =  3/2Z.   However^  the  angle  0   is  now  changed 
and  is  given  by 


(^.21) 


ao 


£n 


^dE. 


idr 


12' 


(l  -  cos  0, 


12 


sin  0, 


1 


evaluated  at  r   ,  r   . 


Setting 


2l+   - 


(4.22) 


dE! 


>^ 


1  2E^(^o-  \^ 


dr 


12 


12 


^120        ^-^ 


*e^        dr 


12 


120 


ve  have 
(4.25) 

and 

(4.24) 


"120 


-  iV^^-^^KV^i^-^J 


[(E   /^>-l] 


o'    1 


cos   G 


20 


r         =  2r 
120  10 


(2   + 

0    ) 

(4  + 

0    ) 

/' 

^0' 

4  +  0 


The 


quantities  E   ,   E     entering  into  0  are  evaluated  at  r  ^    .      Once 


r         is   obtained  by  iteration^   Equation   (4.9)   then  gives 


(1(.25) 


OO  ^  '  \         ' 

^120 


/// 


-2Z(r^+  r^) 


( V  ^2)(1  -    cos  Qg)  2,  ,      ,     ,, 

X     r^  dr^  r^  dPg  sin  ©^  dO^ 

^12 
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The  remaining  integration  can  be  done  in  several  vays.   It  is  simplest  to 
use  the  Hylleraas  variables.   The  integral  then  becomes 


oo 


IT 


dr 


12 


/.J 

^12  ° 


12 


dt   s(r^2-   *)    ^' 


■2Zs 


vhich  is   easily  evaluated  and  equals  5/(2Z)    .      Finally, 


(^.26) 


H 


.A  E    (E   -  Ej 

A  o^    o        1 


oo 


120 


Table  II  presents  a  comparison  of  approximate  and  exact  values  of 
Hqq  for  various  Z.   The  error  is  again  less  than  one  percent  although 
larger  than  in  the  case  of  the  normalization  integral. 


z 

(?)■' 

(<)■- 

00 
00 

(;>j^"^ 

(hzA        'A 

'A 
H     ^ 
00 

^•Ex 

00 

1 

9.^58^ 

8.95I2 

0.9^^2 

5.802^ 

5.^\ 

0.906^ 

2 

5.055^ 

5.0753 

1.0072 

1-370^ 

i.560g 

0.992^ 

5 

2.095-^ 

2.1114 

1.008 

0.977q 

0.980^ 

1.003^ 

k 

1-737^ 

1.7^^6^ 

1.005^ 

0.8265q 

0.8505^ 

1.0052 

5 

1-5555 

1.5586 

1.002^ 
5 

0.7^706 

0.7510 

1.005^ 

6 

l.kk^Q 

l.H5o 

l.OOOg 

0.6982g 

0.7020 

1.005|^ 

7 

1.5696 

1.5693 

0.999g 

0.66580 

0.6690 

l.OOi^^ 

TAJBLE  II. 

A      'A  'W'v  ''Ptt 

Comparison  of  approximate  values.  N  ,  H    vith  exact  values, N   ,  H 

00   00  00    00 


26  - 


Before  closing  this  section,  ve  will  indicate  the  way  in  which  the 
approximation  technique  will  be  applied  in  lithium.   The  function  E^ 
recommends  itself  because  it  is  readily  integrable.   To  obtain  the  best 
accuracy  possible  in  lithium  consistent  with  a  reasonable  amount  of  labor 
we  will  always  integrate  over  one  of  the  interelectron  coordinates  as  well 
as  the  radial  coordinates.   To  see  how  this  is  to  be  done  the  differential 
element  of  volume  is  examined.   Choosing  one  electronic  vector,  say  r  , 
as  the  polar  axis  and  denoting  9  ,  Q  ,  0  as  the  azimuthal  angles  of 
r  ,  r   with  respect  to  r   and  the  polar  angle  of  r   with  respect  to  a 
plane  through  r^  and  r  ,  we  easily  see  that 

(i+.27)    dY  =  8zt^  T^   dr^  v^    dr^  r^  dr  sin  9^  d9^  sin  0^  dQ^  A.(/)^    . 


We  now  choose  r   as  one  of  the  Integration  variables.   From 


^12 


=  /  r^+  r  -  2r  r  ( cos  9  cos  0  +  sin  9  sin  9  cos  0  )  we  obtain 


2     .       .2.       .     ^^1^^2 


(i^.28)     dV  =  8it  r^  dr^  ^2   ^^2  ^3   ^^5  ^12  "^^12 


sin  02 


In  conformity  with  the  previous  discussions  r  ,  r     ,    r  ,  r^^Q  ^^^ 
chosen  by 

(I+.29)  [  f  X  (r  -  r   )  dV  =  0      1  =  1,2,5 


J 


f  X  (r^2-  ^120)  dV  =  0   . 
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To  choose  the  angular  expansion  points,  we  note  that  there  is  a  singularity 
in  the  integral  in  the  0  direction  at  0  =  0,rt  because  of  the  volume 
element  (4.28).   In  practice  the  0  =  «  value  vhich  places  r  and  r 
farthest  away  from  each  other  gives  the  largest  value  in  the  vicinity  of 
the  singularity.   We  thus  take 

(^.30)  i^20  =  "   • 

0   and  0   are  now  determined  by  the  condition  analagous  to  (4.5),  that  is. 


(^.31)  ^  (^  g)  =  0  ;   ^  an   g)  =  0 

subject  to  Equations  (4.29)  and  (4.50).  As  mentioned  an  iteration  procedure 
is  used  to  solve  both  equations  in  (4.51)  simultaneously. 


5.   Energy  calculation  for  the  ground  state  of  lithiiun 

The  wave  function  for  the  ground  state  of  lithium  follows  from 
Equation  (3-6).   The  two  core  electrons  in  lithium  will  be  in  the  (is) 
hydrogenic  state  and  the  outer  electron  will  be  in  the  (2S)  hydrogenic 
state.   The  interelectron  wave  functions  are  all  of  the  same  form.   The 
determinant  indicated  in  (3-6)  is  then  formed  remembering  that  proper  spin 
functions  must  be  included.   Using  the  symmetry  properties  of  the  Hamilton- 

ri5i 

ian  and  of  the  spin  functions^  -■ ,  the  wave  function  leduces  to  the  following 
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form  for  calcalating  the  energy 


(5.1)    1= 


■Zr. 


-  Zrg   -  2  r 
e      e 


(1 


-Zr   --r^ 
e     e 


(1  -  I  r^) 


X  E^  (r^^)  ^1  ^^23^  ^1  ^'" 


25'   1  ^  51' 


=  I  (1,2,5)  -  $  (1,5,2) 


Exchange  effects  occur  because  of  the  interaction  of  the  first  and  second 
terms  on  the  right. 

The  energy  is  obtained  from  Equations  (5-2)  and  f5'5)'   For  the 
ground  state  of  lithium  the  energy  is 


(5.2) 


E   = 


-9Z  (M       +  M        ) 
n,e       e,e 

8  N 


W 


^     dr     dr     dr^ 


(5.5)     M 


n.e 


=  >    >  M. . .  =  />    y   r  ^  cos  Q. . . 


B^i 


1  +  J 


i   ^   J 


Sr.    Sr.  . 
1        ij 


dV 


M 


.e   =     21       ^       51     V, 
k  4=    Jt    4  m 


km 


Jl 


cos    0, 


s^  $ 


kLkm     ~x  N. 

k  +     je      +  m         '  ^^kjg  ^^km 


dV 
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M    expresses  the  coupling  of  the  interelectron  and  nuclear  coordi- 

nates  and  is  much  larger  than  M    which  expresses  interelectron-inter- 

e,  e 

electron  coupling.   In  the  usual  manner  the  M's  and  N  can  be  split  into 
two  parts^  one  of  which  expresses  the  direct  interactions  while  the  other 
expresses  the  exchange  or  overlap  effects.   Thus 


N     =  N^-   N^ 


T)  V 

i,ij  ijij  ijij 


where 


\;?,km         ^  kjp.km  "  ^  kj?,km 


(5.5)  ^  =  2  j  j  jf^  (1,2,3)   dV 

^   =  2  I  I  [$  (1,2,3)  1  (1,5,2)    dV 

D              (   f   f                          a^l,2,5)}     ^ 
M.    .  .   =  cos   0.    .  .     —  dV 


i,iJ        J   )   )  i.iJ  Sr.    Sr.  . 

E  f  (  (  S2{$(i,2,5)  I  (1,5,2)}   ^ 

M.     .  .    =  cos    0,     .  .      ::  dV 

^^^J        J    )    }  i.iJ  ar     5r 

1       ij 


D        r  r  r  ^^^^  (1.2,5)} 

\   ,lan  =  j    J   J    -^   \i,km     7" " 


E                (    (  (                             ^^{1  (1.2,3)   1  (1,3,2)} 
^,lan  =  J  j  J    -^  \^km    T" -T 


dV 
k£       km 

'(£  u,^,3;  i  u.3,i^;^ 

dV 


Br,    «    Sr, 
kii        km 
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Note  that  the  derivatives  act  explicitly  on  the  coordinates  indicated  in 
the  differentiation.  After  dlfferentiating,the  r. .  must  be  expressed  in 
terms  of  r.  and  r   and  the  integration  performed.   The  integration  is 
done  approximately  in  the  manner  outlined  in  the  previous  section  by  inte- 
grating one  of  the  r. .  variables  along  with  the  remaining  radial  and 
angular  coordinates.   The  Integrals  in  the  M's  and  in  N  contain  polynomials 
in  their  Integrands.   The  approximation  technique  is  applied  term  by  term 
to  the  polynomials  in  order  to  retain  the  greatest  accuracy. 

The  actual  integrations  vlll  be  illustrated  in  several  cases.  When 
■written  in  full  N   has  the  form 


,  -2Z(r  +  r  )   -Zr  2 

(5.6)   W^=Je     1   2^   5  (1  _  2^^  ^  Z_^^2)  g2^^^^^  g2^^ 


23' 


X  E^2  (^^^)  ^f 


The  Interelectron  variable  chosen  for  Integration  connects  the  radial  vari- 
ables vhlch  have  the  same  exponential  dependence.   The  choice  is  r   in  this 
case.   The  Hylleraas  variables  can  now  be  used  to  perform  the  necessary 
Integrations.   Take  the  first  integral  on  the  right  as  an  example.   There 


(5.7)  f  =  e-22"  E^(r^2^ 


^=^l'^-23)'l'^V 


dV  =  8jt^  (s^-  t^)  r   ds  dt  dr   r  ^  dr 


-  51  - 


Now 


■2Z(r^+  r^)   2 


(5.8)     J  e  eJ  (^12' ^V  ^10^  '^'^  =  ° 


-        =  -  i  ^  ^  N 
10     2  N    dZ   00 
00 


where  N   is  evaluated  as  before  and  the  derivative  is  readily  performed 
00 

using  Equations  (l.l)  and  (l.2)  in  the  Appendix.  We  note  r   =  r   . 
For  Z  =  5,  r^Q  =  r^^  =  O.5625O.   r   is  given  simply  by 


(5.9)     J   ^     ^'■5"  ""50^  <iV  =  0  ;   r. 


-Zr,  ^ 

■30  -   I  =  1-°°°°° 


The  formula  for  r  _  is 


^   -2Z(r  +  r  ) 
(5.10)    J   e     ^   ^  E,^  (r,J(r,,-  r,^^)  dV  =  0 


12  ^12''  12   120' 


(2/)  r   /2Z(r^+  ^2^  2  _  2  ,„2  .2, 


^120  =  ^f-"  j   "         "  ^12  ^12  (^  -  t  )  ds  dt  du 
00 


iSzLl     J    ^  f  2   ^2J    +  i  I    IZj    +  -^  I^iZj] 
W^^        I    2Z  V  5     5  Z     4^    '  2     5^    '' 


00 


0.86592  . 
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I'he  angles   are  determined  easily  in  this   case  by  symmetry.      0     =  rt 

as  was  shown  in  the  previous  section.      Thus   r. ^    r^^    r,   are  situated  in 

one  place   and  the   angles  9^ ,    0^   of  r^    and  r^  with  respect  to  r_  will  be 

L         d  L  d  ^ 

equal  at  the  extrema  of  the  integrand,  since  the  integrand  is  unchanged 
with  the  interchange  of  G  ,  Op . 


FIGURE  1. 

Configuration  of  three  electrons  in  a  plane, 
illustrating  various  distances   and  angles. 


Figure  1  displays  the  general  configuration  in  space.     When  0^^  =  9 
and  r,  ^  =  r-„  the  distances  r^^.,    ^o:z/^  a^^  obtained  immediately. 

10  d\j  lldKJ         d^yj 


(5-11) 


2  ^  2 

^130  "  ^250 


"^150  =  ^230  =  1.^261^ 


-  33 


Hence 


,       -2Z(r  +  r   )      -Zr         22? 
(5.12)  j      e  12^        5e^2    (^_^^)   E^2^r23)    E2(r^^) 


dV 


If  2i 

=  4jt  E^    (1.4264)   N 


'  °°  (3Z)5 


The  other  integrals  in  N  are  evaluated  in  the  same  manner.   The  only 
change  vhich  occurs  is  the  change  in  r   which  takes  on  values  (4/z)^ 

(5/z)  in  the  tvo  remaining  integrals  respectively. 

F 
The  overlap  integral  N  has  the  form 


/   -2Zr   -3Z(r  +  r  )/2,  2     v 

(5.13)    W  ^  =  j   e    1  e     2   3   ^1  .  I  (^^^  ^^)  ^  ^  ^^^^j 


XE^2^r^2)\'^-25)^l'^V  ""^ 


The  interelectron  distance  is  chosen  to  be  r  •   The  Hylleraas  variables 
are  again  introduced,  and  the  integrals  dealt  with  individually.   The 
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only  point  to  note  is  that  the  I  must  nov  be  evaluated  at  a  nev  value,  Z  , 
given  by  Z  =  ^  Z  =  2.25.   The  value  r   does  not  change  in  the  different 
integrals  while  r   =  r  ,  ^p^r)   '^°*  '^^^   quantities  r   ,  ^i  on  ^-"^^  again 
equal  and  are  obtained  by  the  symmetry  properties  used  in  N  .   One  point 
of  significance  will  be  noted  here.   The  second  integral  on  the  right 
includes  the  term  (r  +  r  )•   The  greatest  accuracy  is  found  when  each  term 
is  treated  independently  by  splitting  this  second  integral  into  two  integrals. 


5^ 


The  r^_,  r^_  are  no  longer  equal  in  each  integral  and  symmetry  is  then 
20   50 

lacking.   Hovever,  in  calculating  W   the  term  (r  +  r  )  was  treated  as 

one  unit  since  W   is  quite  small  compared  to  N   .   Several  cases  are 

illustrated  below  in  which  terms  like  the  one  described  must  be  split  up. 

The  method  used  in  treating  M    is  now  illiistrated.   Actually 

there  are  two  methods  of  dealing  with  M   .   The  first  follows  the  approxi- 

mation  technique  outlined  for  H   in  Section  ^.   The  second  makes  use  of 
^  00 

integration  by  parts  under  the  integral  sign  to  convert  M    into  integrals 
of  the  type  appearing  in  N.   The  integration  by  parts,  however,  has  the 
effect  of  producing  very  large  cancellation  among  the  resulting  integrals 
so  that  the  answer  sought  is  much  smaller  than  the  integrals  themselves. 
This  second  method  is  thus  highly  inaccurate.   However,  the  comparison 
of  answers  given  by  both  methods  furnishes  a  valuable  internal  check  on 
the  accuracy  of  the  approximation  technique  as  used  in  lithium. 

We  illustrate  the  first  method  as  applied  to  M  .   It  is  convenient 
to  combine  terms  having  the  same  value  of  i,j  and  to  treat  these  together. 
Take  M     +  M      as  an  example.   Using  Equations  (5-5)  aJ^d.  (2.5)  the 
integrals  combine  so  that 


r       -2Z{v  +   r  )-Zr   (r  +  r  ) 
(5.1^)     M  ^^^2  .  M  ^^^^   =  Z  J   e     12    3  ^__^  (1  .  ,,,  ,^^) 

2  2  _ 
x((l  -  Zr,.  !p-)[i_  E/(r^^)]  E^^  (r,,)  E^  (r,^)  dV  . 
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Take  the  first  integral  on  the  right.   The  interelectron  variable  chosen 
for  integration  is,  of  course,  -"^lo'   Then  the  Hylleraas  variables  are 
introduced  and 


(5.15) 


f  =  e 


-2Z(r^.  r^)    ^ 


dr 


12 


■^I't- 


12' 


=  ■^l'('23'  ^/("-Jl! 


The  expansion  points  are  detennined  in  a  manner  similar  to  N   .   For  example, 


(5.16) 


dH 


00 


"10        ' 

2H       dZ 
00 


Thus  r^Q  =  r^Q  =  0. 55626,  r   =  1.00000,  r^^^  =  O.91605.   The  angle  0^0  ^^  ^^ 
vhile  fi-,  „/  G^^  are  given  by  symmetry  so  that  Equation  (5.11)  holds  and 

^130  =  ^250  =  1-5932. 
Finally, 


f       -2Z(r  +  r  )-Zr   (r  +  r  ) 
(5.17)    j   e     123    1   ^  (1  - 

^12 


cos  9^^) 


d?^^'^^12) 


X  E^^r   )  E^2(r   )  dV  =  W  E^^  (1.3932)  H^^ 


21 


°°  (5Z)5 


The  same  procedure  applies  to  the  other  integrals  in  M  -.  o  +  M   o 
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The   other  tenns  in  M       do  not  yield  quite  so  easily  to  the   calcula- 
tion,  but  they  Involve  no  essential  difficulty.     As  an  example,   ve  take 


(■5.18)  M^  +  M^  =  M°  +  M° 

\-y.i.o)  2^23  3,23  1,13  3,15 


2,23  3 


.23  =  ^1 


-2Z(r^+  r2)-Zr^ 


dfc  \^  (^2,) 


25 


^'('•i^'^i^-ji' 


X 


\{v^+  r^)(l  -   cos   9      ) 


( 


2 
Z  2 


1  -   Zr,   +     -f^     r.^  ^ 


^3  "     ^     ^3 


(r^-   Tg  cos  Q^j) 


,,2    2       -I 


>•  dV 


The   calculation  of  the   integral 


(5.19) 


J 


-2Z(r^+  r^)-'Lv^        (r^-   r^   cos   Q^^  )     ^   ^ 


"23 


(¥-3^ 


X 


d „  2,         X 

IT-  "^1   ^^23^ 


23 


\'(-12)\'(^3l)    dV 


will  T^e  illiistrated  and  the  other  tenns  may  be  found  in  an  analogous  manner. 


Here 
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(5.20) 


.2Z(r3_.  r^)   -Zr 


^1  ^^12^  ^3 


g  =G?^ 


(r^-  Pg  cos  9g^) 
"25 


_d 2  /    N 


^lir^^) 


where 


%  = 


(r^-  r^  cos  9^^) 


G  = 


dr 


23 


=1^-23) 


\'(^3l) 


The  expansion  points  r   =  t     ,    r     ,    r    are  obtained  from  Equation  (^.29)^ 

•with  the  above  formula  for  f.  y^   is  integrable  once  the  r   are  obtained. 

Since  the  r„,  and  r^ ,  variables  do  not  enter  into  g  symmetrically,  we 
23      13 

cannot  use  the  results  given  earlier  to  determine  the  angles.   To  obtain 
the  angles  and  thus  r  q,  r  ^,   we  arrange  r  ,  v     ,    r   in  a  plane  as 
before.  Noting  that  0^^_  =  0..^^  +  Q^  ^-  where 

d.^\j  jJ-LU     J-dV 


(5.21) 


cos  0 


"120 


•10 


■20 


120 


2r  r 
10  20 


The  maxim-um  of  Hn   g  with  respect  to  0   can  be  obtained.   This  yields 


sin  0 


(5.22) 


tan  0 


120 


250    7   +  cos  0 


120 


vhere 


58 


(5.23)    7 


r   r 

10  30 


"I30 


r  (^ 


150^ 


(5' ^-250)  -0 


r   r 
20  30 


"230 


cos  9^^^  -  (r^o/r^o) 


vhere  ]^   and  J   are  given  by  Equations  (4.l4)  and  (4.22).   Iteration  is 
used  to  solve  these  two  equations  simultaneously  starting  vith  assumed 
values  of  6   ,  ^12,0'  '^^^   calculations  can  be  set  up  to  proceed  quite 
rapidly  on  a  desk  calculator.   Once  r    and  r    are  determined,  the 
integral  in  (5.19)  is  taken  as 


(i.2k)  c|E^==(r^3„)3j^Ef(,^j^,  (Je'     ^ 


-2Z(r,+r^)   -Zr 


2'  e  "5  E 


^^r^^jdr^dr^r^dr^) 


(r     -  r       cos  G     )  , 

>CI  2n    1      -^ ^ ^       sin  G25  d023 

^250 


(.J 


where 


/ 


-2Z(r,+   r^)      -Zr_  ,  ^      ^ 

12  3  (r  -   r     cos   0      )     ^ 

e  e  r     — ^ ^^  dV 

r 
23 


(5.25)        C  = 


-2Z(r  +  r^)      -Zr^  , 

[e  .        ^  .        ^^20-^30^°^   ^23- 


e  r^ 


dV 


230 


is   the   ratio  of  approximate  to  exact   r^_,^  values  as   Z  ->  00  . 

230 
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The  integrations  involved  in  the  above  are  facilitated  through 
the  relations: 


/   (r. -  r .  cos  9.  .) 
(5.26)    ^  .  .   =  — ^ i^  Sin  9.  .  d9.  . 


iJ    ;        J.  10    iJ 

2 
2  /"j 


r. 

1 


=  2  -  3   -  r^  <  r^ 


3  r  '■j-^'i 


The  integral  in  the  denominator  of  Equation  (5.25)  is  given  immediately 
by  Equation  (5.26) .   The  integral  in  the  numerator  involves  further 
integrations.   It  is  evaluated  by  splitting  the  range  of  integration  of 
r  or  r  and  integrating  in  the  usual  manner.   The  results  are  straight- 
forward although  tedious  and  involve  integrals  of  the  type 


(5.27)  J^(x)    =     j 


00 


n     -2Zx    , 
X     e  dx 


.2Zx     (X-       nx^-1  --^-2 

e  I  ^T^  +  — 


+  n(n  -  1)  =•  +   ...   +  - 

{2Zr  {2Zy  ^2ztl 


The  remainder  of  tne  integrals  in  M  .  .,  +  M  ^  ^,  are  found  in  a 

2,23     5,25 

similar  manner  and  M      can  then  be  evaluated.   An  alternate  expression 

n,e 

for  this  sum  is  obtained  from  Equation  (ll.^)  in  Appendix  II. 


ho 


In  lithium  this  yields 


(5.28)    M^^^^  -j^l^^l2^h'(^2,^h'^^31^ 


I  ^'-    - 


(i  .i  .i 
V  ^1   ^2   ^! 


X 


1  -  Zr, 


^2   2 


5        9  2   2 


-»_ 


-2Z(r  +  r  )   -Zr   _ 
X  e  e   -^  dV 


and  also 


(5.29)    M- 


n,  e 


%'(V=l(-25'"^'-3l' 


_  *-  J 


{ 


X  '^  1  -  z 


(ig+ij)    ^2 


^  F-  ^2  ^5 


l-i^( 


5  -  J  Z  (r2+  r  ) 


+    \   Z^  r^  r^ 


■2Zr^  -5Z(r2+  r^)/2 


dV 


The  disadvantage  with  these  equations  Is  that  there  is  substantial  cancella- 
tion among  the  separate  integrals  which  appear.   The  advantage,  of  course, 
is  that  the  separate  integrals  are  of  the  type  dealt  with  in  the  normaliza- 
tion  integrals  N   and  N   and  their  approximate  integration  is  considerably 
simpler  than  the  integrations  dealt  with  in  M     .  Equations  (5.28)  and 
(5.29)  were  integrated  by  expanding  out  all  the  terms  into  a  single  poly- 
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nomial  eind.  treating  each  term  Individually.  Because  of  the  severe  cancel- 
lation, the  accuracy  lost  by  combining  terms  such  as  (r  +  r  )  into  one 
integral  cannot  be  tolerated.   This  means  that  symmetry  is  lost  in  many 
of  the  integrals  and  the  angles  are  evaluated  by  the  necessary  modification 

of  Equations  (5.22)  and  (5.23).   Otherwise  the  technique  is  the  same  as 

D        F 
used  for  N   and  N 

The  results  found  for  M      using  Equations  (5.28)  and  (5»5)  s^re 
0.01612  and  0.0157^5  respectively.   These  results  furnish  an  internal 
test  of  our  approximation  technique.   The  cancellation  involved  in  Equa- 
tion (5.29)  produces  an  answer  which  is  I5  times  smaller  than  the  largest 
integrals  appearing  there.   This  magnifies  the  errors  in  the  individual 
integrals  by  the  factor,  15?  if  the  errors  are  assumed  uncorrelated. 
Now  evaluating  Equation  (5 •5)  in  the  manner  described  above,  the  answer 
and  the  largest  integral  are  about  equal.   There  is  then  no  magnification 
of  the  errors  in  the  individual  integral.   Table  II  indicates  that  the 
errors  in  the  approximation  method  applied  to  the  normalization  integrals, 
which  appear  in  Equation  (5.28),  and  to  the  perturbation  integrals, 
which  appear  in  Equation  (5-5)  ?  are  comparable  In  the  range  of  Z  for 

which  we  are  interested.   Thus  a  good  measure  of  the  error  in  M 

^  n,e 

as  obtained  by  Equation  (5.5)  is  simply  the  difference  in  the  two  results 
quoted  divided  by  I5..  This  is  an  error  of  less  than  0.2  %  which  compares 

quite  favorably  with  the  results  of  Table  II  for  Z  =  3. 

E 
Because  of  the  good  results  found  using  Equation  (5.28),  M   was 

evaluated  by  Equation  (5.29).   The  small  size  of  M   renders  the  errors 

involved  negligible  for  the  accuracy  we  seek. 


-  k2   - 


The  terms  M     and  M     are  quite  small.   Indeed  for  the  particular 
e,e       e,e 

E  D 

prohlem  at  hand  M     can  be  completely  neglected.   The  terms  in  M 
■^  e^  e  e,e 

are  evaluated  with  sufficient  accuracy  hy  taking  the  cos  ©  j^j^  j^  'ts^™ 
completely  outside  the  integral  sign.   Thus,  for  example;, 

(       -2Z(r  +  rj   -Zr  /dE^      \  /  "^1^  ,    ' 


XE^iv^^)    dV 


-2Z(r  +  r  )  -Zr 
-^  C  |^^(r,_)l  E-(r^,^j  1   e        ^     e       ^ 


&  '^3.0))  ^^^.,c)  J 


■51 

2 

"■■12 


<^12 


where 


2     2     2 
T          —  r   +  r 
12     25    31 
(5.52)    2  cos  e^^ji  =  

^12  ^51 


/ 


-2Z(r^+  Tg)   -Zr   ^ 
cos  9^2,51  -  ^  ^V 


C  = 


/ 


-2Z(r  +  r  )   -Zr   ^ 
s     ^   "^  e   ^  dV 


^5 


The  integral  in  the  numerator  is  evaluated  vith  the  formula 


(5.55)    j  cos  ^i^^^a     -     52«5^^^^ 
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dA  is  the  angular  part  of  dV  and  the -0-  g  are  defined  in  Euqation  (5.26). 
The  subsequent  radial  integrations  must  be  split  properly  aad  use  made  of 
Equation  (5. 27).   In  certain  cases  the  Exponential  Integral'-  -'  arises  and 
the  following  properties  are  used  in  evaluating  the  subsequent  integrals. 

00    . 
e 


(5.5^)  ^   dt  =  -Ei(-x) 


dEi(-x)   ^  £^ 
dx        X 


-Ei(-px)  >  -0.577215  -  j2n  X  -  in  p 

X  ->  o 


The  expansion  points  and  angles  are  treated  in  the  same  manner  as  described 
before . 


6.   Results 

The  energy  calculated  for  the  ground  state   of  lithium  as   described 
in  the  previous  section  is   -7.595  ±     O.OO5  A.U.      The  uncertainty  in  the 
energy  value   is  due  to  the  vairious   approximations   outlined.      It  should  be 
noted  that   a  given  error  in  the  individual  integrations   shows   up  In  the 
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energy  reduced  by  a  factor  of  about  thjree.   The  energy  calculated  here  is 
compai^d  in  Table  III  vith  a  number  of  other  calculations  using  BWT  formed 
from  a  product  of  symmetrized  hydrogenic  vave  functions  but  with  different 
nijmbers  of  variable  parameters.   These  results  can  be  obtained  from  the 
work  of  Wilson '--'-'''-I  . 

The  BWF  vith  no  variable  parameters  are  the  hydrogenic  functions 
vith  Z  =  5.   Screening  is  then  introduced  through  the  variable  parameters. 
With  one  variable  parameter  the  inner  and  outer  orbits  are  screened  in  the 
same  vay  vhile  vith  the  tvo  parameters  they  are  screened  individually  and 
vith  three  parameters  the  node  in  the  outer  orbit  is  adjusted  independent 
of  the  screening.   The  seven  parameter  calculation  of  Huang;,  vhlch  intro- 
duced tvo  of  the  thiree  interelectron  coordinates,  and  the  experimental 
value  I-  -'  are  also  given  for  reference. 


Energy  (A.U. ) 

Number  of  Parameters 

BWF               SWF 

0 

-  7.057 

-  7-595  +  0.005 

1 

-  7-289 

2 

-  7.592 

3 

-  7-^18 

7 

-  7.^57 

Experimental 

-  7.^78 

TABLE  III. 

Compaj'ison  of  varioxis  energy  values  calcu- 
lated for  the  ground  state  of  lithium. 
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The  excellent   result   indicated  in  helium  is  fully  realized  in 
lithium.      As  in  helium  (see  Table   l),    the  SWF  function  is  5  times   closer 
to  the  true   ansver  than  the   zero  order  BWF.      It   is  more  than  tvice   as 
close   as  the   one  parameter  fiinction  and  perhaps   a  little  better  than  the 
two  parameter  function.      The  SWF,    of  course,,    has  no  adjustable  para- 
meters . 

We   can  conclude  that  the  Pluvinage   approach  is   fully  applicable 
to  larger  atcanic  systems.      In  addition^    the  mathematical  approximation 
outlined  can  be  extended  without  great  difficulty  to  these   atomic  sys- 
tems so  that  reasonably  exact   calculations   could  be   carried  out.      This 
is  especially  true   of  quantities  which  might  be   calciilated  from  the  wave 
functions  where  errors  of  the  order  of  0.1%   to  1%    are  tolerable.     We 
also  note  that  the   approximation  technique   csin  be   carried  to  higher 
orders,    if  necessary,    to  obtain  requisite  accuracy.     This,   however,  would 
probably  involve   a  considerable  increase   in  labor. 
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APPEKDIX  I 


1.   The  Integrals,  I 


n 


The  integral  I^  is  given  by  Equation  (2.20).   Successive  I  values 
can  be  generated  by  differentiation.   Thus,  in  viev  of  Equation  (2.19), 


(I.l) 


dyz^ 

dZ 


1 

—  I 
2Z 


n+1 


By  means  of  this  forntula  and  Equations  (l.4)  belov,  we  arrive  at  the 
following  expressions: 
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E     designates  the  modified  Cliff ord-Bessel  function  of  order,    v. 

1     2 
The   argument  in  these   functions  in  the   ahove  equations   is  j-  Z   .     E     has 

ri2i 

the  following  properties^  -• 


J /2i  Ji) 
(1-5)  E  (?;)  = 
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(iv^)' 
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n!( V  +  n) I 
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^    '  ^  v+1     v-1     V 


E   =  E  ^T 
V      v+1 


t 

E       E   ^ 
-/   V        v-1 
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E       E 

V  V 


PlqI 

The  following  recurrence  relation  ^  -^    can  also  be  used  to  obtain  I^i 


(1.5)     8z5  I^^2  =  (l2nZ^  +  Q^)  ^n+1  "  l^^^^^'^'^  +  kn  -  2\    I^ 


+  n(n-l)(n-2)  I^_^ 


^9 


By  using  the  recurrence  relation  higher  I  can  be  generated  if  only  the 
three  lowest  I  are  known.   This  can  be  done  quite  rapidly  on  a  desk 
calculator.   However^  accuracy  is  lost  in  evaluating  the  higher  I  by  the 
recurrence  relation  and  there  is  also  a  possibility  of  carrying  along 
errors  in  the  lower  I  .   Values  of  I  used  in  this  paper  were  always 
computed  by  the  explicit  formulas  (l.2)  and  the  recurrence  relation  (I.5) 
and  found  identical  within  the  accuracy  desired  (four  significant  figures), 
This  seirved  as  an  independent  check  on  the  correctness  of  the  formulas  in 
(1.2). 


APPEHDIX  II 


The  alternate  expression  for  M     is  obtained  by  noting 

2        2  2  N 

(II. 1)        cos  0,    -  -  -^ ^ iJ-  -  -^ 


^i  ij  i 


The  derivative  on  the   right  is  the  partial  of  r       with  respect  to  r     keep- 
lag  r.    constant.     Then 


(II.2) 


2  2  2 

1     \  ^i   -   ^j   ^  ^ij  ^     _     ^ 


J^i  ^i   ^ij  ^ij        ^1 


where  the  derivation  on  the  right  applies  to  functionfl  which  depend 
explicitly  on  the  r.      alone.     We  can  then  convert   certain  derivatives  with 
respect  to  the  Interelectron  coordinates  to  derivatives  with  respect  to  the 
radial  coordinates. 
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From  Equations   (5.6)^    (5- 3)    and.  the   above 


(11-5)    M 


n.e 
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IT 
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/  -^%i,  4:  ^'^^^^'^ 


^     \  p  >  q 
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IT 


J-  (°1T//)   ^(tT    -^-p,)) 
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dV 


The  term  0  (r^)  has  been  abbreviated  to  0£  and  D  Indicates  the  Slater 
determinant  of  the  functions  immediately  following  this  symbol.   The 
second  equation  in  (II. 5)  can  nov  be  integrated  by  parts^  by  explicitly 
expressing  r   In  terms  of  r  ,  r  .   The  boundary  terms  drop  out  in  the 

pq         p  q 

integration  by  parts,  leaving 
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J  J  J  J 
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STT0^ 
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The  argument  in  the  E^s  has  been  suppressed  for  convenience.     Prom 
Equation   (ll.4)_,    the  Schroedlnger  Equation: 


1         1 


and  the  form  for  the  vave   f\:nctlon  given  by  Equation   (5.1)^    the  expres- 
sions  (5.28)    and  (5.29)   follow. 


*     *     * 
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